Anomalous fermion bunching in density-density correlation 
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We consider theoretically density- density correlation of identical Fermi system by including the 
finite resolution of a detector and delta-function term omitted in the ordinary method. We find an 
anomalous fermion bunching effect, which is a quantum effect having no classical analogue. This 
anomalous fermion bunching is studied for ultracold Fermi gases released from a three-dimensional 
optical lattices, ft is found that this anomalous fermion bunching is supported by a recent experiment 
(T. Rom et al Nature 444, 733 (2006)). 
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Demonstrated firstly by Hanbury Brown and Twiss 
[H, quantum noise correlation is a fundamental prob- 
lem in identical Bose and Fermi system, and has im- 
portant application in astronomy, quantum optics, con- 
densed matter physics and subatomic physics 0, H, 0, @ 
etc. In the last ten years, quantum noise correlation was 
studied experimentally for cold Bose and Fermi gases 
d, 0, i, II ES EL El EE E| with thc remarkable de- 
velopment of cold atom physics. Together with these 
experimental advances, intensive theoretical studies con- 
tribute largely to our understanding of correlation effect 
for cold atomic system. The theoretical studies about 
high-order correlation have extended from a single har- 
monically trapped Bose gas EH to different matter state 
such as pair correlations of fermionic superfluid 17 1 
and different trapping potential such as high-order corre- 
lation of cold atoms in an optical lattice [18[ . The finite- 
temperature effect [3| and low-dimensional effect 
about high-order correlation were also studied theoreti- 
cally. 

The theoretical studies for fermionic superfluid lil 17 1 
and one-dimensional ultracold Fermi gas [2l| have shown 
that quantum noise correlation provides important in- 
formation about the fundamental quantum features of 
ultracold Fermi gases. Most recently, there is a clear ob- 
servation of fermion antibunching in a degenerate Fermi 
gas released from a three-dimensional optical lattice [3] ■ 
The fermion antibunching physically arises from the an- 
ticommutation relation of field operators, and thus has 
no classical analogue. It is a manifestation of Pauli ex- 
clusion principle in high-order correlation. Besides the 
observation of the antibunching effect in density-density 
correlation for fermionic atoms released from an optical 
lattice El, fer mion antibunching was also observed for 
electrons 22 , 23 , 24 1 and neutrons [2|| . 

In the present Letter, we find that when the resolu- 
tion of a detector for the measurement of density dis- 



tribution is considered, there would be an anomalous 
fermion bunching effect under appropriate conditions. 
This anomalous fermion bunching effect originates phys- 
ically from the delta function in the anticommutation 
relation of fermion field operators, omitted in the ordi- 
nary method in calculating the density-density correla- 
tion. The consideration of finite resolution of the detector 
will avoid the notorious divergence in the delta-function 
term, and lead to anomalous bunching effect. The ob- 
servation of obvious anomalous fermion bunching effect 
requests special conditions and parameters. Fortunately, 
we notice that there is already an anomalous fermion 
bunching in the experimental data of Ref. EH ■ 

Because we will give a comparison between our the- 
ory and experimental data to support the anomalous 
fermion bunching, we first introduce briefly the elegant 
experiment in Ref. [Hi]. In this experiment, an ultracold 
Fermi gas of K was firstly prepared in the combined po- 
tential of an optical trap and a 3D (three-dimensional) 
optical lattice. After free expansion of 10 ms by sud- 
denly switching off the combined potential, 2D (two- 
dimensional) density images were recorded with a CCD 
(charge-coupled device) camera by illuminating a reso- 
nant laser along the vertical (z) direction. See Fig. 1(a). 
The 2D density-density correlation was obtained by deal- 
ing with appropriately a set of 2D density images. 

The starting point of our theory is the 2D density- 
density correlation function given by 
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(1) 

Here n is a 2D density operator, r = {x, y} and 
d = {d xi d y }. To deal with appropriately the layered 
distribution of fermionic atoms in 3D optical lattice and 
absorption imaging along z direction, the 2D density op- 
erator takes the form n(r,t) — gW (r , t)^ (r , t) . The op- 
erator g has the property that ^ = f(dij)- 
gij = J^k fijkz w hh fijk z being the occupation number 
in the lattice site indexed by {i, j, k z }. The summation is 
about the lattice site along z direction, g^ represents the 
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overall particle number in a string of lattice sites along 
z direction, is an annihilation operator for the atom 
in an equivalent 2D lattice site indexed by {i,j}- The 
2D density-density correlation function is dependent on 
d. C*2(d,t) < 1 corresponds to the fermion antibunch- 
ing effect, while C<i(d, t) > 1 corresponds an anomalous 
fermion bunching effect. 

We give here a brief introduction about the reason for 
the above rules in the transformation from 3D system to 
2D density-density correlation measured by a CCD. Be- 
cause the density images were recorded in Ref. [l3} with 
a CCD camera by illuminating a resonant laser along z 
direction, special consideration should be given about the 
2D density-density correlation for a 3D system. The 2D 
density distribution is 
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Here <pij (r, t) is the wave function of the atom initially 
in the lattice site Because in 2D density-density 

correlation, the coordinate z has been integrated, it is 
convenient to define the 2D density operator as 



n(r,t) = <7$ f (r,i)$(r,i). 



(3) 



With the above rules about g, (n(r, t)) is the same as the 
result given by equation ((2]). With this definition of 2D 
density operator and <?, it is straightforward to get the 
final result of the 2D density-density correlation function. 
From Eq. fl}, it is easy to get 
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The delta-function term in the above formula is or- 
dinarily omitted because the divergent prop- 
erty and <5(d) = for d ^ 0. This term originates 
from the delta function in the anticommutation rela- 
tion |$(r + d/2,i),$ t (r - d/2,*)} = 5(d), and thus ac- 
counts for a pure quantum effect. 

Bunching (antibunching) corresponds to a peak (dip) 
in density-density correlation. For d — > 0, the delta- 
function term in equation (j4|) means a divergent bunch- 
ing behavior. When the finite width Ad of spatial resolu- 
tion is considered, it is understandable that there would 
be an effect of increasing the width and decreasing the 
height (depth) in the peak (dip) of the bunching (an- 
tibunching) behavior. Roughly speaking, 6(d) may be 
replaced by a function with height 1/A? in the region 
-Ad/2 < d x < A d /2 and -Ad/2 < d y < A d /2, and 



with zero value outside this region. When both the delta- 
function term and spatial resolution are considered, more 
accurate results are obtained by calculating 

«n(r- d/2, t)n(r + d/2, *)»,,= 



J d 2 s 1 J d 2 s 2 (n(r - d/2 + Sl , f)n(r + d/2 + s 2 , t)} . 

(5) 

Here S denotes the region —Ad/2 < s x < Ad/2 and 
—Ad/2 < s y < Ad/2. In the above formula, ()d repre- 
sents the average due to the finite resolution of the de- 
tector. This average has been considered in the starting 
point Q. 

From equation |T]), we get the following approximate 
result 
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Here I = 2Ttht/ml p with m being the atomic mass and 
l p being the spatial period of the optical lattice. S is 
the area of the overlapping region between Si (deter- 
mined by -Ad/2 < x < A d /2 and -A d /2 < y < A d /2) 
and S2 (determined by —Ad/2 < x — d. x < A d /2 and 
—Ad/2 < y — d y < Ad/2). N is the total particle num- 
ber. In equation ((6]), the second term represents the 
antibunching behavior. This term is obtained by omit- 
ting the effect of spatial resolution and under the con- 
dition l/a > Lo/lp and I ^> A d l p /a. Here L and a 
represent respectively the overall width of the Fermi gas 
and wavepacket width of an atom in a lattice site before 
switching off the combined potential. The last term in 
equation ^ physically originates from the delta- function 
term in the anticommutation relation of field operators. 
This term is obtained by considering the spatial resolu- 
tion and under the condition l/a > Lq/I p . 

Using the experimental parameters in Ref. [13j and 
equation ©, Fig. lb gives C2(d,i) for flight time of 10 
ms. Besides eight dark dots representing antibunching, 
at the center location we notice a bright dot representing 
bunching. Fig.Ic gives further one of the dark dot, while 
Fig. Id gives further the bright dot. Eight dark dots are 
due to the second term in equation @. This term is —I 
for d x /l = i and d y /l = j with i and j being integers. At 
the locations of these dark dots, the last term in equa- 
tion ([6|) is zero. The bright dot at the center is due to 
the last term in equation (|6|) , which is larger than 1 close 
to d x — and d y = 0. If the last term in equation ([6]) 
is not included, there should be a dark dot at the center 
(see Fig.le). In Ref. the theoretical model with- 

out the last term is used to interpret their experimental 
data. It is quite interesting to note that eight dark dots 
(rather than nine dark dots without the consideration of 
the last term in equation (O) were observed in Ref. [l3| 
(see Fig. 2c and its figure caption in this reference)! In 
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FIG. 1: Fermion antibunching and bunching in density- 
density correlation function, (a) The measurement of density 
distribution for a Fermi gas released from a three-dimensional 
optical lattice, (b) Two-dimensional density-density correla- 
tion obtained from equation ((6} and experimental parameters 
in Ref. [l3j]. The theoretical results agree with the exper- 
imental data. In particular, an anomalous fermion bunch- 
ing is shown clearly by the bright dot at the center, which 
has been supported strongly in Ref. [l^ . (c) Illustration of 
a dark dot (antibunching effect) in Fig. lb. (d) Illustration 
of the bright dot (bunching effect), (e) The density-density 
correlation near the center without the consideration of the 
delta-function term in the ordinary theory. In these figures, 
the spatial coordinate is in unit of fim. 

Fig. 2d of Ref. [13;], the obvious and "mysterious" peak 
at the center of density-density correlation shows further 
the anomalous fermion bunching behavior. This com- 
parison gives strong experimental evidence for anoma- 
lous fermion bunching. It is clear that the condition for 
observing fermion bunching is l 2 l 2 T,{jg 2 j > 2-nb? d a 2 N 2 . 

The measurement of density-density correlation ex- 
tracts the correlation information in quantum noise. We 
discuss further the physical mechanism of the bunching 
effect for Fermi system through the consideration of den- 
sity fluctuations. The density fluctuations are defined 
as 

5 2 n(x, t) = lim (ra(x, t)n{y, t)) - (n(x, i)) 2 . (7) 

y^x 

Simple calculations give 

«S 2 ?i(x, t) = lim (*t( x t )$( Vj 4 )\ s( x - y) - (n(x, t)) 2 . 

y^x \ / 

(8) 

We see that there is a divergent <5-function term in 
5 2 n(x, t). It is obvious that the omission of the 5-function 
term will lead to absurd result of negative density fluctu- 
ations. The divergent (^-function term is a pure quantum 
effect by noting that it comes from the anticommutation 
relation between field operators. There is another origin 
for this (^-function term that the field operators ^(x, t) 
and tyi(x,t) comprise infinite modes, rather than only 
the modes occupied by particles before a measurement. 
In fact, this property is an essential property of quantum 



field theory. The analyses about 8 2 n(x, t) show clearly 
that the J-function term can not be omitted simply. Sim- 
ilarly to our studies of the density-density correlation, the 
divergence in the J-function term can be avoided because 
the resolution of a detector always has a width. This is 
equivalent to the presentation that creation or annihila- 
tion of a particle at an infinitesimal point is impossible 
because this would mean an infinite energy exchange. 

Note that the anomalous fermion bunching does not vi- 
olate in any sense the Pauli exclusion principle for identi- 
cal fermions. The anomalous fermion bunching originates 
from the delta-function term in the anticommutation re- 
lation of field operators. Of course, the Pauli exclusion 
principle will try to destroy the anomalous bunching ef- 
fect. The second term in ^ reflects the Pauli exclusion 
principle, and leads to the antibunching effect, and even 
can completely destroy the bunching effect at the center. 

In summary, our studies give an anomalous fermion 
bunching, and we have given a strong evidence by ana- 
lyzing a recent experimental data in Ref. [13;]. The con- 
dition to observe the anomalous fermion bunching is in 
fact quite rigorous. It is not surprising that this unique 
quantum effect is found accidentally in Ref. [3] with 
highly developed experimental technique. The experi- 
mental technique in Ref. jl3[ gives us chance to test 
further our theory, (i) The last term in equation ([6]) 
increases with the increasing of flight time. Thus, we ex- 
pect a transition process from the antibunching to bunch- 
ing at the center of C2(d, t), with the increasing of the 
flight time, (ii) The distribution of fermionic atoms in the 
lattice sites {i,j, k z } is determined by i 2 +j 2 + a 2 z k 2 < R 2 
at zero temperature. Here a z is determined by the har- 
monic potential due to the optical trap. Simple calcu- 
lations give Hygfj/N 2 ~ 0.93(a z N)- 2 / 3 . We see that 
decreasing a z has an effect of enhancing the anomalous 
fermion bunching effect for identical particle number. 
This dependence on a z would give us further chance to 
test our theory. We believe further experimental stud- 
ies would deep largely our understanding of the indis- 
tinguishability of identical particles, wave packet local- 
ization in quantum measurement process and high-order 
correlation. High resolution of a detector is required to 
reveal anomalous fermion bunching effect. The astonish- 
ing resolution of STM (scanning tunneling microscope) 
and AFM (atomic force microscope) etc suggests that 
the experimental and theoretical studies about electrons 
in an atom or molecule may lead to a new regime about 
the studies of high-order correlation, quantum noise and 
quantum measurement etc. 

The ^-function term in the anticommutation relation 
is the direct reason for the anomalous fermion bunching 
effect. For Bose system, the inclusion of the 5-function 
term in the commutation relation of field operators may 
also play important role in density-density correlation. It 
is easy to understand that the inclusion of the 5-function 
term will lead to an enhanced boson bunching effect at 
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d = 0. In fact, in a recent experiment about the density- 
density correlation of ultracold bosonic atoms released 
from an optical lattice, there is a clear enhanced bo- 
son bunching effect at d = in C 2 (d, t) (see Fig. 2(c) 
in Q, where the brightness in the center dot is much 
higher than other eight bright dots, and this can not 
be explained without the consideration of the ^-function 
term.). This sort of enhanced boson bunching effect is 
also implied strongly in [12j. These analyses show that 
anomalous fermion bunching effect and enhanced boson 
bunching effect have common physical origin — the indis- 
tinguishability of identical particles and the 5-function 
term in the anticommutation (commutation) relation of 
field operators. 
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